Effects of Quark Spin Flip on the Collins Fragmentation Function in a Toy Model 
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The recent extension of the NJL-jet model to hadronization of transversely polarized quarks al- 
lowed the study of the Collins fragmentation function. Both favored and unfavored Collins fragmen- 
tation functions were generated, the latter purely by multiple hadron emissions, with 1/2 moments 
of opposite sign in the region of the light-cone momentum fraction z accessible in current exper- 
iments. Hints of such behavior has been seen in the measurements in several experiments. Also, 
in the transverse momentum dependent (TMD) hadron emission probabilities, modulations of up 
to fourth order in sine of the polar angle were observed, while the Collins effect describes just the 
linear modulations. A crucial part of the extended model was the calculation of the quark spin flip 
probability after each hadron emission in the jet. Here we study the effects of this probability on 
the resulting unfavored and favored Collins functions by setting it as a constant and use a toy model 
for the elementary single hadron emission probabilities. The results of the Monte Carlo simulations 
showed that preferential quark spin flip in the elementary hadron emission is needed to generate 
the favored and unfavored Collins functions with opposite sign 1/2 moments. For the TMD hadron 
emission modulations, we showed that the model quark spin flip probabilities are a partial source of 
the higher rode modulations, while the other source is the Collins modulation of the remnant quark 
from the hadron emission recoil. 

PACS numbers: 13.60.Hb, 13.60.Le, 13.87.Fh, 12.39.Ki 
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I. INTRODUCTION 

The recent interest in the transverse momentum struc- 
ture of hadrons has prompted renewed efforts both in ex- 
perimental and theoretical studies. Notably, experiments 
using semi-inclusive deep inelastic scattering (SIDIS) 
have advanced our understanding of the underlying par- 
tonic structure of the nucleon in momentum space. For 
example, the naively time-reversal odd Sivers distribu- 
tion function has been measured to be non-zero [TH3] . In 
the factorization regime, the SIDIS cross-sections can be 
expressed as convolutions of parton distribution functions 
with elementary parton-probe scattering amplitudes cal- 
culable using perturbative QCD and parton fragmenta- 
tion functions ^4; .5, ■ The complete tree- level expressions 
at leading order have been presented in Ref [5]. Thus, 
to access the underling partonic structure of the nucle- 
ons, a precise knowledge of the fragmentation functions 
are needed. The particularly interesting case of the T- 
odd Collins fragmentation function has been widely 
explored both in theoretical models [7HTD] and experi- 
ments [U El [TTHT^ , with the first direct measurements of 
the Collins fragmentation mechanism being performed 
by the Belle collaboration using hadron pair production 
in e"^e~ collisions [Tl] [T^]. Further, the experimental 
studies from HERMES, COMPASS and JLab are sug- 
gesting that the unfavored Collins functions have a sim- 
ilar size and an opposite sign to that of the favored 



ones [DiaiTllTl]. 

The theoretical calculation of Collins function for pi- 
ons and later also kaons within the spectator model of 
Refs. [THSI made predictions for the favored Collins func- 
tion using the mechanism of interference of one-meson- 
loop amplitudes with the tree level amplitude. However, 
these model calculations could not make direct calcula- 
tions of the unfavored Collins function, as only a single 
hadron emission was included. Recently, the Nambu- 
Jona-Lasinio (NJL) jet model of Refs. [T5HT^ has been 
extended to calculate the Collins fragmentation func- 
tions, both favored and unfavored, using Monte Carlo 
(MC) simulations in the quark-jet hadronization mecha- 
nism [20]. Here the transverse polarization of the frag- 
menting quark has been introduced, and the light cone 
spinors were used to calculate the quark spin flip proba- 
bilities in each hadron emission step. In the hadroniza- 
tion of a transversely polarized quark, the Collins effect 
describes a modulation of the unpolarized hadron frag- 
mentation function with a Collins fragmentation function 
term that is proportional to the sine of the angle ip of the 
hadron's transverse momentum and the quark's spin in 
the 7*A^ frame, as depicted in Fig. [l] The correspond- 
ing expression using the "Trento Convention" [3T], with 
polarized quark q carrying momentum k and spin Sq 
fragmenting to unpolarized hadron h of mass m/j, carry- 
ing the light-cone momentum fraction z and transverse 
momentum P± with respect to quark's momentum k can 
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FIG. 1. Illustration of the three dimensional kinematics of 
transversely polarized quark fragmentation. The fragmenting 
quark's momentum fc defines the z-eods with its transverse 
polarization spin vector Sq along x axis. The emitted hadron 
has momentum p with the transverse component Pj_ with 
respect to the z-axis. The polar angle of hadron's momentum 
P with respect to the zx plane is denoted by ip. 



be written as 
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where the unpolarized fragmentation function is denoted 
dI'''{z,PI) and H^'^'^z.Pl) is the Colhns function. 
When integrated over P|^, the polarized fragmentation 
function can be expressed in terms of the integrated un- 



polarized fragmentation function {z) and the 1/2 
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In the Monte Carlo simulations of Ref. |20], we used 
as input the elementary unpolarized and Collins frag- 
mentation functions for single hadron production, calcu- 
lated within NJL model, where we used the mechanism 
of Ref. |F for the elementary Collins function. The full 
polarized fragmentation function was then calculated in 
MC simulations with multiple hadron emission. The re- 
sulting polarized fragmentations, when integrated over 
the hadron transverse momentum squared (Pj), exhib- 
ited the Collins modulation both for favored and unfa- 
vored hadron production. Here the unfavored Collins 



function's 1/2 moment is comparable in size but oppo- 
site in sign to that for the favored fragmentation. For 
the transverse momentum dependent (TMD) polarized 
fragmentation functions, at fixed values for the hadron's 
light cone momentum fraction z and , exhibited mod- 
ulation with a fourth order polynomial in sin((/3), even in 
simulations with only two produced hadrons. 

In this work we have two goals, focusing on pion frag- 
mentations from light quarks only. First, we examine the 
dependence of the sign and the magnitude of the 1/2 mo- 
ment of unfavored fragmentation function on the quark 
spin flip probability used in the NJL-jet model MC sim- 
ulations. Second, we examine the origins of the higher 
order Collins modulations in the TMD polarized frag- 
mentation functions. 



II. NJL-JET MODEL WITH TRANSVERSELY 
POLARIZED QUARK 




FIG. 2. NJL-jet model including transverse momentum and 
quark polarization transfer. Here the orange double-lined ar- 
rows schematically indicate the spin direction of the quark in 
the decay chain. 

Schematically, the NJL-jet model with the transversely 
polarized quark is depicted in Fig. [2] In MC simula- 
tions of Ref. [2D], in each hadronization step, we used 
the elementary polarized fragmentations to sample the 
momentum of the production hadron, zi, p± and ipi, 
with respect to the fragmenting quark in that step (here 
ipi is defined with respect to the spin of the fragmenting 
quark). Using the momentum we calculated and recorded 
the hadron's momentum, z, P± and (p, with respect to 
the initial fragmenting quark. Further, we also deter- 
mined the momentum of the remnant quark using mo- 
mentum conservation, while the probability of the rem- 
nant quark's spin flip was sampled using the spin non-flip 
and flip probabilities |aip and |a_ip. These probabilities 
were calculated using Dirac spinors of the transversely 
polarized quark, that in turn can be expressed in terms 
of the Lepage-Brodsky spinors in helicity basis [22l [23] . 
Explicit these probabilities are proportional to 

|ai|2 ^ ll |a_i|2 ^ ll + (M2 - (1 - z)Mi)2, (5) 

where Mi and M2 are the masses of the fragmenting 
and the remnant quarks, and l± = (Ix^ly) = —p± is the 
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transverse momentum of the remnant quark with respect 
to the momentum direction of the fragmenting quark. 

In this work, to explore the role of the spin-flip proba- 
bility, we perform MC simulations with various fixed val- 
ues of this probability. This allows us to clearly demon- 
strate the role of the quark spin flip during the hadroniza- 
tion process on the resulting fragmentation functions. 
Also, using constant for this probability instead of the 
one with the dynamical factors of Eq. ([s]) will help to 
disentangle the sources of the higher order modulation 
effects in sin((/3) observed in TMD fragmentations. 

For input to the Monte Carlo framework, the elemen- 
tary probability of a transversely polarized quark q to 
emit a single hadron h carrying light-cone momentum 
fraction z, transverse momentum p± (with magnitude p± 
and polar angle ip with respect to the emitting quark's 
spin), can be written as 

-i?i^''/^(z,pi)^sinM, (6) 
zruh 

where d^/'^{z,p\) is the elementary unpolarized frag- 
mentation function and H^'^^'' {z,p\) is the elementary 
Collins fragmentation function. The unpolarized frag- 
mentation function has been calculated within the NJL 
model in tree level approximation [18|. The correspond- 
ing expression for an unpolarized quark emitting a pseu- 
doscalar meson can be written as 

d\'''{z,p\) = GhQ{pi)j^glqQ z (7) 

pj -f [{z - l)Mi + Ma]^ 
[p]_ + z{z - l)Ml + zM| + z)mlf ' 

where the subscripts on the constituent masses, Mi and 
M2, denote the flavors of the fragmenting and the rem- 
nant quarks respectively, that are also indicated by the 
subscripts q and Q. and ghqQ denote the isospin 
factor and quark- meson coupling constant, that have 
been determined within the NJL model [HI [TB], while 
GhQ{p\) is the dipole regulator for the transverse mo- 
mentum dependance. The details of the derivations and 
the model parameters can be found in Ref. [18] 

The model elementary Collins function within the NJL 
model has been also calculated in Ref. [20]) using the 
mechanism of interference between tree level amplitude 
and amplitude with gauge link coupling for single hadron 
emission, outlined in Refs. [HI l^- Here we are inter- 
ested in exploring the general properties of the quark-jet 
hadronization mechanism, thus the details of the par- 
ticular models are unimportant. To further simplify the 
model and help reduce the required number of the Monte 
Carlo events that allow for a reliable extraction of the 
higher order modulations by enhancing the azimuthal 
modulation signal with respect to the statistical fluctua- 
tions, we employ a toy model for the elementary Collins 
function taking it promotional to the unpolarized one. 



with a coefficient chosen such that the polarized proba- 
bility is still positive for all values of the corresponding 
arguments 

= - 0.9sin^), (8) 

where we assume that H^^^^ (z.p^, )^J^ = 

0.9 di^'^{z,p]_) and Sq = 1. Within the NJL-jet model, 
the elementary splitting functions are renormalized such 
that quark's total probability of emitting a hadron in 
each step is one: / dz dp\/2 dtp di^/q^{z,p']_, (p) = 1, 
the sum is over all the hadrons, while a quark of given 
fiavor can emit directly. These renormalized splittings 
will be used in the next two sections as input to the 
Monte Carlo simulations of the quark-jet hadronization 
process. The renormalized elementary unpolarized 
fragmentation function, integrated over p^, for frag- 
mentation of a M quark to 7r+ is depicted in Fig. [s] All 
the other elementary fragmentation functions from light 

quarks (u and d) to pions are related to d^ either by 
a simple isospin factor or isospin symmetry. 
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FIG. 3. Elementary renormalized unpolarized fragmentation 
function, used as input in the Monte Carlo simulations. 



III. THE QUARK SPIN FLIP PROBABILITY 
EFFECTS ON 1/2 MOMENTS OF THE COLLINS 
FRAGMENTATION FUNCTIONS 

In this section we explore the properties of the 1/2 mo- 
ments of the Collins functions, Eq. Q, for pions emit- 
ted in the hadronization of an up quark as a function of 
the quark spin fiip probability SF in each hadron emis- 
sion step. We perform Monte Carlo simulations to ex- 
tract the polarized quark fragmentation functions D/i/gt 
as functions of z, and tp within the NJL-jet frame- 
work, as described in Ref. [20]. In this section, we limit 
the number of emitted hadrons in each quark decay chain 
to Nj^inks = 2. This allows us to cleanly see the effects of 
the quark spin flip on the emission of the second hadron 
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4. Fitted values for twice the Collins function 1/2 mo- 
(a) and u — >■ 7r~ (b) as a function 



ment, 2H^^^^'^^ for u — >■ tt ' (aj ana u — >■ tt 
of z from Monte Carlo simulations with three different values 
of the remnant quark spin flip probability SF in each elemen- 
tary emission, where the number of emitted hadrons in each 
quark decay chain is fixed to N Links ~ 2. 



in the decay chain, where the unfavored fragmentation 
functions are first generated. 

Next, to extract the integrated unpolarized fragmen- 
tation function ([3| and the 1/2 moment of the Collins 
fragmentation functions Q, we use the relation to the 
P]_ integrated transversely polarized fragmentation writ- 
ten in Eq. ([2|. Thus, we perform a minimum-^^ fit to 
Df^jq\{z,ip) using a form F{cq,ci) = cq + cisai{tp) for 
fixed values of z, identifying the fitted coefficients with 
the unpolarized and Collins terms in Eq. ^ . The result- 
ing fits for the integrated pion polarized fragmenta- 
tion functions with the functional form of F(co,ci) per- 
fectly describe the results of the MC simulations, with 
average per number of degrees of freedom (NDF) very 
close to unity. The unpolarized fragmentations are not 
affected by the quark spin fiip probability SF and have 



been studied in detail in our previous work |171ll8j . Thus 
we only present here the results for the Collins functions. 

The results from fits for H^'"'^^'^\z) for fragmentations 
of a u quark to tt"*" and ■n" are depicted in Fig. |4j where 
each curve depicts the results for a different value of 
SF. The results for 7r+ in the subplot (a) show that, if 

SF < 0.5, then ff^^^^^^ is always positive, crossing zero 
and becoming negative at small values of z for SF > 0.5. 
For unfavored tt~ fragmentation, the effects of SF are 
more dramatic: H^'^^^'^^ is positive (same sign as that for 
7r+) for SF < 0.5 and vanishes for SF = 0.5. Thus, in 
order to generate a non-zero unfavored Collins fragmen- 
tation function, the quark spin non-flip and flip must be 
unequally probable. Finally, for 5*^" > 0.5 the H^^^^'^'^ 
is negative, at 5*^" — 1 becoming exactly the negative of 
the results for SF = 0. This is easy to understand, as 
the unfavored functions are only generated after the ini- 
tial hadron emission, so we can consider only two hadron 
emission for simplicity. In case where SF = 1 (flip of the 
quark's spin after each emission), the angle in Eq. ([6]) 
that is measured with respect to the spin of the initial 
quark acquires an additional phase tt compared to the 
simulations with SF — (no quark spin flip) since the 
modulation in each hadron emission step is generated 
with respect to the azimuthal angle of the currently frag- 
menting quark, resulting in the Collins term having the 
opposite sign. 

Thus, the representative features of the full model re- 
sults in Ref. [ID], where the 1/2 moments behave similarly 
to the scenario where SF > 0.5, qualitatively do not de- 
pend on the details of the particular input unpolarized 
and Collins fragmentation functions, but rather represent 
the characteristics of the quark-jet model with _preferen- 
tial spin flip probability (|a_ip > 



\aiY in Eq. 5 



IV. HIGHER ORDER MODULATIONS IN TMD 
COLLINS FUNCTION 



In the previous section we showed that setting quark 
spin flip probability SF = 1 can be considered as the 
asymptotic solution for generating the largest value of 
the unfavored Collins function that mimics the behavior 
of our model calculations in Ref. [20] . Thus in this section 
we will concentrate on simulations with SF — 1, that will 
allow us to easily extract the Collins modulations from 
the unpolarized fragmentation functions. 

In this section we consider the full transverse momen- 
tum dependent fragmentation functions. As mentioned 
in Ref. |20| . the polarized fragmentation function can no 
longer be successfully described with a first order poly- 
nomial in sin((p) for all values of z and Pj, even for MC 
results with just two hadrons emitted, Nunks = 2. The 
plots in Fig. [5] depict the MC result D^+/y_t{z,P'^,ip) 
(in blue dots) for a fixed value of z = 0.3 and = 
0.04 GeV^ (a), and P| = 0.25 GeV^ (b). The minimum- 
fits with a second order polynomial in sm{(p) are de- 
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(a) 




(b) 



5. The histograms (blue dots) for the polarized number 



FIG. 

density -D^+/„t {z, PI,v) for fixed values of Pi = 0.04 GeV^ 
(a) and = 0.25 GeV^ (b) with z — 0.3 for both cases, as 
a function of the azimuthal angle ip, from Monte Carlo sim- 
ulations using the toy model with a fixed number of hadrons 
emitted in each decay chain, NLmks ~ 2. The minimum-x^ 
fits with a second order polynomial in sin((/p) (red lines) are 
also depicted. 



picted as red lines. 

Here we analyze the dependence of the order of these 
modulations on the number of hadrons produced, N Links- 
For that, we consider minimum-x^ fits of the correspond- 
ing polarized fragmentation functions with polynomials 
in sin(</7) of various orders, for fixed values of z and 
Pj. To quantify the quality of the fits, we use as a 
measure the /NDF, that should be close to 1 for a 
good fit. Our MC simulations are performed with 100 
discretization points for both z and in the interval 
[0, 1) and 100 discretization points for if in the interval 
[0, 27r). Thus, the number of the "data points" to fit, for 
all distinct values of z and P]_, as well as all the possi- 
ble pion fragmentations of u quark (to 3 types of pions) 
are Np = 3 x 10'*. The histograms in Fig. 6] depict the 
results for these fits for simulations with Nunks =2 (a), 
N Links = 3 (b) and N Links =6 (c). Here the horizon- 
tal axis denotes the /NDF and the vertical axis is the 



number of the fits that yield this value, while the differ- 
ent lines represent the results for fits with different order 
polynomials in sin((^). The results for N Links = 2, de- 
picted in Fig. |6ja) , show that the linear form of Eq. ([T]) is 
inadequate to describe the polarized fragmentations for 
a significant number of fixed values of z and P^. The 
quadratic form, on the other hand, provides an excellent 
description of the simulations, and it can be confirmed 
by the histograms that cubic fits don't improve the re- 
sulting values of chi"^ /NDF. Similarly, for the simulation 
with N Links = 3 of Fig. [6]^b) , the cubic polynomial yields 
excellent fits, while the lower order polynomials are in- 
adequate to fully describe the simulations. Finally, for 
the simulations with N Links = 6, depicted in Fig. [6]jc), 
at least a fourth order polynomial is needed for good fits, 
though there is a marginal improvement with fits using 
sixth order polynomials. 

A simple explanation of these modulations within the 
quark-jet picture can be easily given. Here, after the 
initial hadron emission, the hadron emission probabil- 
ity is modulated according to the elementary polarized 
fragmentation of Eq. ([S]). Thus, the remnant quark, hav- 
ing transverse recoil momentum with respect to the ini- 
tial quark that is equal to minus that of the emitted 
hadron, also acquires opposite signed sine modulation in 
the distribution of its transverse momentum. In the next 
hadron emission step, the emitted hadron probability is 
again modulated according to the elementary polarized 
fragmentation of Eq. (|3|, but now with respect to the 
currently emitting quark, that itself has a sine modula- 
tion in its polar angle distribution. Thus, the produced 
hadron probability acquires a second order sine modula- 
tion in if with respect to the initial quark, as does the 
remnant quark in this step, and so on. Thus, in the 
simulations with Nunks produced hadrons, the resulting 
polarized fragmentation functions should be modulated 
with a polynomial of order N Links , as we have seen from 
the results of Figs. [6] (a) and (b), while the simulations 
with Nunks — 1 yield the trivial linear modulation. We 
have to also note from our earlier discussions of the mech- 
anism for this quark recoil traverse momentum induced 
modulations, that the hadrons emitted in n-th step will 
only have a emission probability modulation of the order 
n in sin((/j). Also, we have shown in our previous work, 
that the remnant quark looses almost all of its light-cone 
momenta within just a few hadron emission steps |17) . so 
that the hadrons produced with higher and higher orders 
of modulations will have smaller and smaller average val- 
ues of z. This explains the fact that the linear form gives 
reasonably good fits for significant number of slices even 
in the simulations with large number of N Links- More- 
over, the results in Fig. for Nunks = 6 also can be 
explained using this argument, as our simulations, with 
only 100 uniform discretization points in z, do not probe 
in detail the very low z region, where the sixth order 
polynomial modulations will be significant. Thus the cor- 
responding results can be well described with just fourth 
order polynomials. Also, the size of the statistical errors 
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FIG. 6. The histograms of the values Nf of xV^DF for 
all the hadrons produced by a w quark with polynomials of 
different orders in sin(iy9) for MC simulations for values of 
N Links of 2 (a), 3 (b) and 6 (c). 



and relatively small number of the discretization points 
in if would affect the sensitivity to this higher order terms 
in the fits. 

In the full model calculations, where the quark spin flip 
probabilities were taken to be those of Eq. ^ , fourth or- 
der modulations of the polarized fragmentation functions 
appear in simulations with just two produced hadrons in 
each quark decay chain. The additional two orders of 
modulation in sui(ip) originate from the quark spin flip 
probabilities themselves, as those depend on the squares 
of the components of the recoil quark transverse momen- 
tum, that can be easily expressed in terms of sin^ {ip) 
using the transverse momentum conservation relation 
of the produced hadron p± and the remnant quark 
with respect to the fragmenting quark: l± = —p± and 
Px = Picos^((p), Py = p^sin^lip). In practice, we saw 
that the modulations of order higher than fourth are very 
hard to extract from our simulations with the present 
discretization of variables and the number of MC events 
calculated (of order 10^^ in current work and 10^^ in 
Ref. [20]). This is a purely computational obstacle, that 
can be easily solved by running the simulation on an or- 
der of magnitude more computational nodes than used 
for this work (only 50 Intel Core-i7 cores ran for about 
two weeks), although it is hard to justify the need for 
such a computationally expensive exercise. 



V. CONCLUSIONS 

In this article we examined, using a toy model, the 
role of the quark spin flip in the quark-jet hadronization 
picture on the pion Collins functions produced in Monte 
Carlo simulations. This work was put forward to com- 
pliment and elucidate the findings for the 1/2 moment 
of the Collins function in our previous work in Ref. [20] 
and also to provide the toy model analyses needed for 
explanations of the higher order sine modulations in the 
transverse momentum dependence of the polarized frag- 
mentation functions, first observed in the same work. 

First, the dependence of the 1/2 moments of the 
Collins functions, H^^'^^^\ on the quark spin flip prob- 
ability (SF) for simulations with only two hadron emis- 
sions was examined in Section IIIII It was shown that 
^_jJ-(i/2) j^g^g ^ very strong dependence on SF, both for 
favored and unfavored fragmentation functions. In par- 
ticular, a non-zero unfavored i/j^^^^^' is only possible if 
SF ^ 0.5, with the results for SF = 0.5 ± a, for any 
< a < 0.5, having equal magnitude and opposite sign. 
The results with SF = 1 have the largest magnitude 
and opposite sign to the input elementary favored Collins 
function, similar to the results of Ref. [5D]. For the fa- 
vored fragmentation functions, the results for SF = 
are of the same sign as the elementary favored Collins 
function for all values of z. As the SF increases from 
to 0.5, the low-z region gets increasingly suppressed. For 
SF > 0.5, they start positive (with respect to the ele- 
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mentary Collins function), then cross zero and become 
negative at small values of z, with minimum negative 
value decreasing as SF increases towards 1. The results 
for SF > 0.5 again resemble those acquired in the full 
model calculations of Ref. [20 , dictating that the qual- 
itative features of H^'^^^'^^ are generated by the quark 
spin flip process in the hadronization process and don't 
depend on the particular elementary Collins functions 
used as input to the MC simulations. 

In Section |IV| we examined the higher order modula- 
tions of the transverse momentum dependent polarized 
quark fragmentation functions for SF = 1, in order to 
determine the role of the model spin flip probabilities of 
Eq. ([5]), used in Ref. [20j, in the mechanism that gener- 
ates these modulations. We studied the order of the poly- 
nomials in sin((^) needed to describe the polarized frag- 
mentation functions from MC simulations with various 
numbers of the produced hadrons Nunks- We concluded 
that the order of these modulations is equal to Nunks 



when SF is kept fixed, though the numerical examina- 
tions of modulation orders larger than 4 were limited by 
the variable discretization and statistical fluctuation ef- 
fects. This part of the modulation was attributed to the 
modulation of the remnant quark in the decay with the 
recoil to the hadron emissions. Thus, the modulations 
that rise by additional two orders with a unit increase in 
N Links observed in Ref. [20j are attributed to the form of 
the model SF that depends on the components of trans- 
verse momentum of the remnant quarks. 
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